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Abstract 

We present exact solutions for the zero-temperature partition function (chromatic polyno- 
mial P) and the ground state degeneracy per site W (= exponent of the ground-state entropy) 
for the g-state Potts antiferromagnet on strips of the square lattice of width Ly vertices and 
arbitrarily great length L^. vertices. The specific solutions are for (a) Ly = 4, {FBCy, PBCx) 
(cyclic); (b) Ly = 4, {FBCy.TPBC^) (M6bius); (c) Ly = 5,6, {PBCy.FBC^) (cylindrical); 
and (d) Ly = 5, {FBCy,FBC^) (open), where FBC, PBC, and TPBC denote free, peri- 
odic, and twisted periodic boundary conditions, respectively. In the L^ oo limit of each 
strip we discuss the analytic structure of W in the complex q plane. The respective W 
functions are evaluated numerically for various values of q. Several inferences are presented 
for the chromatic polynomials and analytic structure of W for lattice strips with arbitrarily 
great Ly. The absence of a nonpathological L^ ^ oo limit for real nonintegral q in the 
interval < g < 3 (0 < g < 4) for strips of the square (triangular) lattice is discussed. 
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1 Introduction 



The g-state Potts antiferromagnet (AF) [|l], |^ exhibits nonzero ground state entropy, 5*0 > 
(without frustration) for sufficiently large g on a given lattice A or, more generally, on a 
graph G. This is equivalent to a ground state degeneracy per site W > 1, since 5*0 = 
fcelnVr. Such nonzero ground state entropy is important as an exception to the third law 
of thermodynamics 0. There is a close connection with graph theory here, since the zero- 
temperature partition function of the above-mentioned g-state Potts antiferromagnet on a 
graph G satisfies 

Z{G,q,T = Q)pAF = P{G,q) (1.1) 

where P(G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color 
(for reviews, see 0-[§l). The minimum number of colors necessary for such a coloring of G 
is called the chromatic number, xiG). Thus 

W{{G},q)= hm P(G,g)^/" (1.2) 

where n = v{G) is the number of vertices of G and {G} = lim„^oo G. At certain special 
points qs (typically = 0, 1, .., x{G)), one has the noncommutativity of limits 

lim lim P(G,g)^/" ^ lim lim P(G,g)i/" (1.3) 

q— >(js n— >oo n— >oo (J— x/s 

and hence it is necessary to specify the order of the limits in the definition of H^dG}, g^) [0- 
Denoting Wgn and Wng as the functions defined by the different order of limits on the left 
and right-hand sides of ( p.. 3D , we take W = Wgn here; this has the advantage of removing 
certain isolated discontinuities that are present in Wng- 

Using the expression for P{G, g), one can generalize g from Z+ to C. The zeros of P{G, q) 
in the complex g plane are called chromatic zeros; a subset of these may form an accumulation 
set in the — > oo limit, denoted B, which is the continuous locus of points where ({G}, g) 
is nonanalytic. [| The maximal region in the complex g plane to which one can analytically 
continue the function W{{G},q) from physical values where there is nonzero ground state 
entropy is denoted Ri. The maximal value of g where E intersects the (positive) real axis is 
labelled gc({G}). This point is important since it separates the interval g > gc({G}) on the 
positive real g axis where the Potts model (with g extended from Z_|_ to M) exhibits nonzero 
ground state entropy (which increases with g, asymptotically approaching Sq = /c^lng for 
^For some families of graphs B may be null, and W may also be nonanalytic at certain discrete points. 
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large g, and which for a regular lattice A can be calculated approximately via large-g series 
expansions) from the interval < g < qc{{G}) in which 5*0 has a different analytic form. 

In the present work we report exact solutions for chromatic polynomials P{G, q) for strips 
of the square lattice with arbitrarily great length L^. vertices of the following types: (a) width 
Ly = A vertices and {FBCy, PBC^) = cyclic; (b) Ly = A and (FBCy, TPBC^) = M5bius, 
(c) Ly = 5 and {PBCy, FBC^) = cylindrical; and (d) Ly = 5 and {FBCy, FBC^) = open, 
where FBC, PBC, and TPBC denote free, periodic, and twisted periodic (i.e. periodic 
with reversed orientation) boundary conditions, respectively. For each of these, taking the 
infinite-length limit, we calculate the degeneracy per site, W{{G},q), and the continuous 
nonanalytic locus B. A comparative discussion is given of these results together with previous 
exact solutions for strips of smaller widths. These strips of regular lattices are examples of 
recursive families of graphs, where the latter are constructed by successive additions of 
subgraph units to an initial subgraph. 

There are several motivations for this work. We have mentioned the basic importance of 
nonzero ground state entropy in statistical mechanics [^. Physical examples are provided 
by ice [§-[10| and certain other hydrogen-bonded molecular crystals |jri|. From the point 
of view of rigorous statistical mechanics, exact solutions are always valuable for the insight 
that they give into the behavior of the given system under study. Although infinite-length 
finite-width strips are quasi-one- dimensional systems and hence (for finite-range spin-spin 
interactions) do not have finite-temperature phase transitions, their zero-temperature critical 
points are of interest. Indeed, the presence of a zero-temperature critical point for the Ising 
antiferromagnet0 on infinite-length, finite- width strips of the square lattice has an interesting 
connection with the behavior of the singular locus B for the strips that we have studied with 
global circuit^: in these cases, this singular locus passes through the point g = 2 in the 
complex q plane. Our exact solutions for P and, in the ^ oo limit, W , thus show 
quantitatively the relation between critical behavior as a function of temperature (at T = 0) 
in the free energy and singularities as a function of q in the per-site ground state degeneracy 
W . The present results also show many interesting connections with mathematical graph 
theory, as is clear from the identity (|1 . 1|) , and algebraic geometry, as follows from the fact 
that for these strips, B is an algebraic curve. Besides the works already cited, some related 

^We recall that on bipartite graphs such as cyclic strips of the square lattice with even Lx, an elementary 
mapping shows the Ising ferromagnet and antiferromagnet to be equivalent; since the Lx ^ oo limit can be 
taken with even Lx, this implies that the critical behavior of the Ising ferromagnet is equivalent to that of 
the Ising antiferromagnet on infinite-length limits cyclic strips of the square lattice. By similar elementary 
reasoning, one can show this equivalence for the infinite-length limit of Mobius strips of the square lattice. 

global circuit is a route following a lattice direction which has the topology of the circle, , and a 
length £g.c. that goes to infinity as n — *■ oo. For strip graphs, global circuits are equivalent to periodic or 
twisted periodic boundary conditions. 
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work on chromatic polynomials for recursive graphs includes |jT2|-||3^; further discussion of 
background and references may be found in [Q. 

A generic form for chromatic polynomials for a strip graph of type Gs, width Ly, and 
length Lj. is 

Nx 

P{G,,Ly X L,,BCy,BC,,q) = ^ c,(g)(A,(g))^^ (1.4) 

i=i 

where Cj{q) and the A^a terms Xj{q) depend on the type of strip graph Gg including the 
boundary conditions but are independent of L^. 

2 Ly = A Square-Lattice Strips with {FBCy, (T)PBC^) 

In this section we give our solutions for the chromatic polynomials of the Ly x strips 
of the square lattice with {FBGy, PBGx) and {FBGy,TPBGx), i.e. cyclic and Mobius, 
boundary conditions, respectively. For both the cyclic and Mobius strips, for > 3 to 
avoid certain degenerate cases, the square lattice strips of width Ly have n = L^Ly vertices 
and e = L^i^Ly — 1) edges. The cyclic square strips have x = 2 for L^ even and x = 3 for 
Lx odd, independent of Ly. For Mobius square strips, x = 2 for {Lx,Ly) = (e, o) or (o, e), 
and X = 3 for (L^., Ly) = (e, e) or (o, o), where e and o denote even and odd. 

We calculate the chromatic polynomials by iterated use of the deletion-contraction theo- 
rem 0, together with coloring matrix methods [|l^, The calculation is considerably more 



involved than that for the Ly = 3 cyclic strip given in p8[ , as is indicated by the number of 



Xj terms in eq. ( |1.4|) , namely, A^a = 26, as contrasted with the value A^a = 10 for the Ly = 3 



cyclic strip. Elsewhere we have given a general determination of A^a as a function of Ly 



As Ly increases, the number of terms A^a in (|1.4| ) grows rapidly; it is 70, 192, and 534 for 



Ly = 5, 6, and 7. We obtain the exact solutions of the form (|L 

26 

P(sg(4 X m, FBGy, PBG,), q) = J2 c.g4,,(A.,4,,)'" (2.1) 

and 

26 

P{sq{4 X m, FBGy, TPBG,), q) =J2 Csg4M^i(A.,4j)'" (2.2) 

i=i 

where L^ = m. The fact that the Aj's for a Mobius strip must be the same as those for the 
cyclic strip of the same width and lattice type was proved in |2^; this also proves, a fortiori. 



that (i) the total number, Nx, of Aj's, and (ii) the continuous nonanalytic locus B, including 
the point qc, are the same for the cyclic and Mobius strips of a given type. For B and qc. 
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we shall often indicate this by the notation {FBCy, {T)PBCx)- The explicit A^g^'s that we 
calculate are as follows. The first six are 

Asg4,l = 1 (2-3) 

Asg4,2 = 3 - g (2.4) 

Asg4,3 = 1 - g (2.5) 

Asg4,(4,5) = (3 ± ) - g (2.6) 

and 

A.,4,6 = (g-i)(g-3) . (2.7) 

Three of the remaining Asg4j's, labelled j = 7,8,9, including the one that is dominant in 
region Ri, are identical to the three that enter into the chromatic polynomial for the Ly = 4 



strip with {FBCy, FBCx) boundary conditions, which was previously calculated in 
This identity was shown in [Q. The remaining Asg4j's for 10 < j < 26 are roots of another 
cubic equation, for j = 10, 11, 12; a quartic equation for 13 < j < 16; and two 5th degree 
equations, for 17 < j < 26. Since the equations defining these A^-'s are somewhat lengthy, 
we give them in the appendix. In Table |l] we list various properties of our calculation and 
compare them with the properties that we have found for other related strips of the square 
(and triangular) lattice. The results for the Ly —>■ oo limit for the triangular lattice with 



{PBCy, FBCx) are from |T^. Comparisons for other lattices such as honeycomb and kagome 



were given in |21|, |23|, |2|, 

In particular, the fact that for this width, Ly = 4 for the cyclic strip of the square lattice, 
we encounter equations of degree 5 for the Aj's means that it is not possible to solve for the 
corresponding A/s as algebraic roots. Our experience with lattice strips of a given width Ly 
(and arbitrary length) and a given set of boundary conditions is that the maximal degrees 
of the factors in the general equation for the Aj's are non- decreasing functions of Ly. Thus, 
assuming that this property of non-decreasing degrees of algebraic factors in the equation 
for the Aj's continues for higher Ly, our present results indicate that the exact solutions in 
| 28| for the Aj's for the width Ly = 3 strip of the square lattice have completed the program 



of obtaining exact algebraic expressions for these terms for this type of lattice strip. 

Although no closed- form algebraic expression can be obtained for the Aj's, a theorem 
on symmetric polynomials of roots of algebraic equations, discussed in [^, enables one 



to calculate the chromatic polynomials exactly to arbitrary order. The key to this is the 
property that since the chromatic polynomial for a cyclic strip is a symmetric polynomial in 
the various roots, it can be expressed in terms of the coefficients of the algebraic equations 
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Table 1: Properties of P, W, and B for strip graphs Gg of the square (sq) and triangular (tri) lattices. New 

results in this work are marked with an asterisk in the first column. The properties apply for a given strip of 
type Gg of size Ly x L^, some apply for arbitrary L^, such as Nx, while others apply for the infinite- length 
limit, such as the properties of the locus B. For the boundary conditions in the y and x directions {BCy, 
BCx), F, P, and T denote free, periodic, and orientation-reversed (twisted) periodic, and the notation (T)P 
means that the results apply for either periodic or orientation-reversed periodic. The column denoted eqs. 
describes the numbers and degrees of the algebraic equations giving the \G,,j\ for example, {3(1), 2(2), 1(3)} 
indicates that there arc 3 linear equations, 2 quadratic equations and one cubic equation. The column 
denoted BCR lists the points at which B crosses the real q axis; the largest of these is for the given 
family Gg. The notation "none" in this column indicates that B does not cross the real q axis. The notation 
"int;(ji; (7c" refers to cases where B contains a real interval, there is a crossing at qi, and the right-hand 
endpoint of the interval is qc- Column labelled "SN" refers to whether B has support for negative Re{q), 
indicated as yes (y) or no (n). 
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that determine these A/s. However, the fact that it is no longer possible to calculate the X/s 
as algebraic roots when the width of the cyclic square strip is 4 means that the determination 
of the nonanalytic locus B must be done in a somewhat more cumbersome manner than in 
our previous work where we had exact algebraic expressions for these Aj's. 

The coefficients Cj that enter into the expressions for the chromatic polynomial ( p..4|) for 
the cyclic and Mobius strip of the square lattice of width Ly are certain polynomials that 
we denote c^'^\ given by [p^] 



k=l 

where 



^j^j for k = l,2,...d (2.9) 

with < d < Ly. We hst below the specific c'^'^^'s that appear in our results for the Ly = A 
square lattice strip: 

= 1 , c(^) = g - 1 , c(2^ = - 3g + 1 , (2.10) 

c(3) =g3_5^2^g^_^ ^ (2.11) 

and 

= (g-l)(g3_6g2 + 9^_l) _ (2.12) 

In ascending order of degrees of c^'^\ we calculate 

Csqi,j = c^°^ for 6 < J < 9 (2.13) 

c^^^ . = c(i) for 13 < j < 21 (2.14) 

c^g^ . = c(2) for 10 < j < 12 and 22 < j < 26 (2.15) 

c^g^ . = c(3) for 2 < J < 5 (2.16) 

and 

Cs,4,i = c(^) . (2.17) 



In it was shown that the coefficient for the Xj that is leading in region Ri must be 1. 
We define 



C{G) = Y,CG,. (2.18) 



6 



where the G-dependence in the coefficients is indicated exphcitly. Note that for recursive 
graphs hke the strip graphs considered here, the depend on Ly and the boundary con- 
ditions, but not on L^. Our results above give 

C{G) = q{q - \f for G = sq{Ly = 4, FBCy, PBC^) . (2.19) 

in accord with the generahzation [03, 123 



C{Gs{Ly X L,,FBCy,PBC,),q) = P(T,.^,g) = g(g - 1)^--^ (2.20) 

for Gs a strip of the square (or triangular) lattice, where P{Tn, q) is the chromatic polynomial 
for the tree graph T„. This is in accord with the coloring matrix approach jS^, |30[| . 
For the Ly = 4 Mobius strip of the square lattice, we find 

Cs,4Mb,j = c(°) for 7 < J < 12 (2.21) 

CsqAMhj = -c^°^ for j = 6 and 22 < j < 26 (2.22) 

CsqAMbj = c^^^ for 17 < j < 21 (2.23) 

CsqAMbj = -c^^^ for J = 1 and 13 < j < 16 (2.24) 

Csq4Mb,j = c^^) for j = 4,5 (2.25) 

and 

CsqAMb,j = -c^^^ for j = 2, 3 . (2.26) 
Hence, the sum of the coefficients is 

C{G) = for G = sq{Ly = 4, FBCy, TPBC^) (2.27) 

in accord with the general result for the Mobius strip of the square (and triangular) lattice 
3| 



^ j P{TLy+i,q) for odd 

> CG{Ly,Mb),j = i n ' f T ■ ^^■'^^> 

^-^ 10 for even L,, 

Chromatic zeros for the cyclic strip of the square lattice with = 4, L^^ = m = 20 and 
hence n = 80 are shown in Fig. |l|; with this value of m, the complex chromatic zeros lie 
close to the boundary B and give an approximate indication of its position. Note that there 
is a zero very close to g = 2, but P{sq{Ly x L^, FBGy, PBCx),q) is nonzero for g = 2 for 
the case shown, where = m is even, as is clear from the fact that x = 2 in this case. 
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Figure 1: Chromatic zeros for the Ly = 4, = m = 20 (i.e., n 



cychc strip of the square lattice. 



The maximal point at which B crosses the real axis, Qc, is determined as a solution of the 



equation of degeneracy of leading terms |A 



eq7—9,max \ 



lA 



eq22—2(3,max \ i 



where A 



eq7—9,max 



and 



Kq22-26,max are the roots of eqs. ( |8.1.1| ) and ( |8.1.5| ) with the largest magnitudes, respectively. 
Since only two Aj's are degenerate in magnitude at this point, it is a regular point on the 
algebraic curve B in the terminology of algebraic geometry. This is also the case for the 
Ly = 3 (and Ly = 1) strip of the square lattice |^8|, |2^, whereas, in contrast, qc is a multiple 
point on i3 for = 2. We find 



qc ~ 2.4928456 for sg(4 x oo, FBCy, {T)PBC^ 



(2.29) 



This may be compared with the values qc = 2 for the Ly x oo strip of the square lattice 
with Ly = 1,2 and the same {FBCy, {T)PBCx) boundary conditions [[7|, and the value 
qc ~ 2.33654 for = 3 |2§. We calculate that W{sq,4:X oo, FBCy, BC^) = 1.2697336.. at 
the value q = qc iia eq. ( |2.29| ). 

The locus B also crosses the real q axis at g = 2 and at g = 0. In addition to region 
Ri which extends outward from the envelope of B and includes the real axis for q > qc and 
g < 0, there are two other regions that contain segments of the real axis: R2, including the 
interval 2 < q < qc and -R3, including the interval < g < 2. In region Ri, the dominant Xj 
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is the root of the cubic equation (|8.1.1|) with the largest magnitude, which we label Xr-g^max- 
In region R2, the dominant Xj is the root of the fifth-degree equation (|8.1.5| ) with the largest 



magnitude, which we label X22~26,max- In region R^, the dominant Xj is the root of the 
fifth-degree equation ( ^.1.41) with the largest magnitude, which we label Xij_2i,max- We have 

= (A7-9,ma.)'/' , for g G i?i (2.30) 
\W\ = \X22-2e,max\'^\ for q G R2 (2.31) 

\W\ = \Xir.2i,max\'^\ for qeRs (2.32) 

(In regions other than Ri, only the magnitude \W\ can be determined unambiguously 0.) 

The locus B has support for Re{q) < as well as Re{q) > 0. It separates the q plane into 
several regions, including the three described above and two complex-conjugate ones which 
we denote R4 and Rl, centered approximately at g ~ 2.6 ± 0.8i. In the regions R4 and Rl, 
we have 

= |Ai3-i6,ma.r/' , for qeR^.Rl (2.33) 

Just as complex-conjugate pairs of tiny sliver regions were found for the cyclic Ly = ?> square 
28| and triangular |^ strips, so also these may be present here; we have not carried out a 



search for such regions (but have ruled out the possibility of tiny regions on the real axis). 

3 = 5, 6 Square- Lattice Strips with (PBCy^ FBC^ 

Here we report our exact solutions for the chromatic polynomials for the width L^, = 5, 6 
strips of the square lattice of arbitrary length and with [PBCy, FBCx), i.e., cylindrical, 
boundary conditions. Results for the cases Ly = 3 and Ly = 4 were given previously in 
H, H, m. We recall that A^a = 1 for Lj, = 3 and Nx = 2 for Ly = 4. For Ly = 5 and 
Ly = 6 we calculate Nx = 2 and A^a = 5, respectively. As before, it is convenient to present 
the results in terms of a generating function, denoted T{Gs, q, x). The chromatic polynomial 
P{{Gs)mi<l) is determined as the coefficient in a Taylor series expansion of this generating 
function in an auxiliary variable x about x = 0: 

00 

T{Gs, q,x) = J2 PiiGs)m, q)x'^ ■ (3.1) 

m=0 

The generating function T{Gs, q, x) is a rational function of the form 
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with 



and 



j=0 



X-' 



V{G,,q,x) = l + Y,bG.A<iy 

i=i 

where the Ac^^i and fec^.i are polynomials in q, and dj^ = degx{J\f), dv 
factorized form 

dx> 

ViGs,q,x) = llil-XG^,iq)x) . 
Equivalently, the Ag^j are roots of the equation 

e^V{G,, q, 1/0 = e- + J2 ^G.d''''' ■ 
The general formula expressing P{Gm, q) in terms of these quantities is 



(3.3) 

(3.4) 
degr,(V), In 

(3.5) 



p{Gm,q) = Yl E^^^: 



dAf 



^d-p—s—l 



j=l ^s=0 



n 



^^i^dx>,ij^j 



For L„ = 5 we find 



A 



1 r 



sq5PF,j 



(A, -A.) J 



J = 1,2 



A-. 



where 



and 



Tsq5PF ± a/ Rsq5PF 

Tsq5PF = q^ - lOg^ + 46g^ - 124g^ + 198g - 148 
Rsq5PF = q^° - 20g^ + ISSg^ - 1092/ + 4356g^ - 12596g^ + 27196g^ 



-44212g^ + 52708g^ - 41760g + 16456 



(3.6) 



(3.7) 



(3.8) 
(3.9) 



(3.10) 

The coefficients Cgq^ppj can be computed using eq. ( p.7| ) in terms of the generating function, 
which is given in the appendix. 

In the La; ^ oo limit, the locus E includes five arcs, consisting of two complex-conjugate 
pairs and a fifth, self-conjugate, arc. The endpoints of these arcs are located at the five 
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Figure 2: Locus B for the width Ly = 5 strip (tube) of the square lattice with {PBCy,FBCx) boundary 
conditions. Thus, the cross sections of the tube form pentagons. For comparison, chromatic zeros calculated 
for the strip length = m + 2 = 1& (i.e., n = 80 vertices) are shown. 
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Figure 3: Locus B for the width Ly = 6 strip (tube) of the square lattice with {PBCy,FBCx) boundary 
conditions. Thus, the cross sections of the tube form hexagons. For comparison, chromatic zeros calculated 
for the strip length = m + 2 = 1& (i.e., n = 96 vertices) are shown. 
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complex-conjugate pairs of roots of Rsqbpp- The self-conjugate arc crosses the real axis at 
the real zero of T^q^pp, namely at 

2.691684 for sq{5 x oo, PBCy, FBC^) (3.11) 

In Fig. ^ we show a plot of chromatic zeros for the Ly = 6 strip of the square lattice with 
[PBCy, FBCx) and length L^. = m + 2 = 16 vertices, so that the strip has n = 96 vertices 
in all. With this large a value of m, the complex chromatic zeros lie close to the boundary B 
and give an approximate indication of its position. (We have not searched for very minute 
features in B.) From our exact analytic results, we calculate indication of its position. From 
our exact analytic results, we calculate 

g^^ 2.6089 for sq{6 x oo, PBCy, FBC^) (3.12) 

The morphology of chromatic zeros for this long 6 x 16 cylindrical strip is similar to that 
found for a 8 x 8 patch of the square lattice, again with cylindrical boundary conditions, in 
P^ . In both cases, the chromatic zeros have support for Re{q) < and prongs extending 
to the right; further, our exact calculation shows that in the limit L^. oo with Ly = 6, 
the locus B has support for Re{q) < 0. In Fig. ^, one of the chromatic zeros is very close 
to g = 2, but for q = 2 exactly, the chromatic polynomial is nonzero, equal to 2, in accord 
with the fact that this strip is bipartite for any value of L^. 
For the ^ oo limit of these respective strips we have 

Wisqib X oo, PBCy, FBC^), q) = (A,,5PF,i,™ax)'/' (3.13) 

and 

W{sq{& X oo, PBCy, FBC,), q) = (A,,6PF,i,™ax)'/' (3.14) 

where \sq5PF,j,max and XsqQPF,j,max dcuotc the solutions to the respective equations ( p.6|) with 
maximal magnitude in region Ri. 

It is of interest to use this exact result to study further the approach of W to the limit for 
the full infinite 2D square lattice. This extends our previous study in |^|. In Table |^ we list 
various values of W{sq{Ly x oo, PBCy, BC^), q) (which, for this range of q, are independent 
of BCx), denoted as W{sq{Ly), P, q), together with Monte Carlo measurements of W for the 
full 2D square lattice, W{sq,q) from ^ and the g = 3 value W^(sg, 3) = (4/3)^/2 f^om [|. 
We also list the ratio 

fl„,(A(L,).i?C,.,)^ ^'^<^^^]^J'-^' (3.15) 

for the present square lattice A = sq. One sees that for Ly = 5 and moderate values of 
g, say 5 or 6, the agreement of W{sq{Ly),q) for the infinite-length, finite-width strips with 
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Table 2: Values of W{sq{Ly x oo, PBCy, BCx,q) (for any BCx), denoted W{sq{Ly), P,q) for short, with 
W{sq,q) — W{sq{oo x oo),q) for 3 < q < 10. For each value of q, the quantities in the upper line are 
identified at the top and the quantities in the lower line are the values of Rwisq{Ly), PBCy, q). 



Q 


W{sq{3),P,q) 


W{sq{4),P,q) 


W{sq{5),P,q) 


W{sq{6),P,q) 


W{sq,q) 


3 


1.25992 
0.8183 


1.58882 
1.032 


1.43097 
0.9294 


1.56168 
1.014 


1.53960.. 
1 


4 


2.22398 
0.9516 


2.37276 
1.015 


2.31865 
0.9921 


2.34339 
1.0027 


2.3370(7) 
1 


5 


3.17480 
0.9766 


3.26878 
1.0055 


3.24518 
0.9982 


3.25196 
1.0003 


3.2510(10) 
1 


6 


4.14082 
0.9858 


4.21082 
1.002505 


4.19790 
0.9994 


4.20058 
1.0001 


4.2003(12) 
1 


7 


5.11723 
0.9904 


5.17377 
1.0013 


5.16557 
0.9997 


5.16689 
1.0000 


5.1669(15) 
1 


8 


6.10017 
0.9930 


6.14792 
1.0008 


6.14221 
0.9999 


6.14296 
1.0000 


6.1431(20) 
1 


9 


7.08734 
0.9947 


7.12881 
1.0005 


7.12458 
0.9999 


7.12506 
1.0000 


7.1254(22) 
1 


10 


8.07737 
0.9957 


8.11409 
1.0002 


8.11083 
0.9998 


8.1111 
0.9999 


8.1122(25) 
1 



the respective values W{sq, q) for the infinite square lattice is excellent; the differences are 
of order 10~^ to 10~^. As noted before PBI, for PBCy (and any BCx) this approach is not 
monotonic. 



4 Ly = b Square- Lattice Strips with [FBCy^ FBCx) 



We have also gone beyond the previous studies in ^ to calculate the chromatic poly- 
nomial for the strip of the square lattice with width Ly = 5 and {FBCy, FBCx), i-e., open, 
boundary conditions. A related study on wide strips is in In ||2^, a given strip {Cs)m 
was constructed by m successive additions of a subgraph H to an endgraph J; here, I = H, 
so that, following the notation of [^, the total length of the strip graph {Cs)m is = m + 2 
vertices, or equivalently, m + 1 edges in the longitudinal direction. The results are conve- 
niently expressed in terms of the coefficient functions in the generating function, as discussed 
above. 

For the width Ly = 5 strip of the square lattice we find degx(T>) = Nx = 7. The coefficient 
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functions bsq5FF,j in eq. (|3.4| ) that determine the Xsq5FFjs via eq. (|3.6|) are hsted in the 
appendix. Because the Agg^ppj's (cf. eq. (p.3|)) are quite lengthy, we do not give them 
here0. In Table |I| this result is compared with the findings from the previous calculations 
in for narrower open strips of the square lattice, and with strips of the triangular 

lattice [f 
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One observes that the equation ( p.6|) defining the \/s increases in degree 
as Ly increases for the open strips. In particular, because we now encounter an equation 
(^]6D of degree higher than 4 (specifically, degree 7), it is not possible to solve for the A/s 
as algebraic roots. Furthermore, assuming that this increase in degree of ( p.6|) continues 
for greater widths Ly of open strips, our present results show that the previous calculations 
of the Aj's in [23| up to Lj^ = 4 have completed the program of calculating these terms 
exactly as algebraic roots for open strips of the square lattice. As noted above, because of 
the theorem on symmetric polynomial functions of roots an algebraic equations p9i , one can 



still calculate the chromatic polynomial in terms of the coefficients of the algebraic equation 



for the Aj's. 
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Figure 4: Chromatic zeros for the Ly — 5 open strip of the square lattice of length — m + 2 — 16 
vertices (i.e. total number of vertices n — 80). 

In Fig. we show a plot of chromatic zeros for the open strip of the square lattice with 



*The Asq5FF.j are listed in the copy of this paper in the cond-mat archive. 
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Ly = 5 and length = m + 2 = 16 vertices, so that the strip has n = 80 vertices in all. 
With this large a value of m, the complex chromatic zeros lie close to the boundary B and 
give an approximate indication of its position. From an analysis of the degeneracy of leading 
Aj's, we find that (in the oo limit where B is defined) 

2.42843 for sq{5 x oo, FBCy, FBC^) (4.1) 

This is in agreement with the chromatic zeros shown in Fig. ^. Comparing Fig. ^ with the 



corresponding plots for Ly = 2 and Ly = 3 (Fig. 3(a,b) of [^), we see that the arcs forming 
B are elongating and that the arc endpoints nearest to the origin are approaching more closely 
to the origin. This agrees with the behavior that we had observed earlier from narrower strips 



and with the conclusions that were drawn from that behavior |2^, |2^, ^ , in particular, 
the statement that these results are consistent with, and provide further support for, the 
hypothesis that in the limit as Ly —>■ oo, the locus B will extend all the way through the 
origin of the q plane and will separate this plane into different regions containing the real 
axis. 

For q > qc, we have, for the physical ground state degeneracy per site of the g-state Potts 
antiferromagnet , 

W{sq{5 X oo, FBCy, FBC,), q) = {\sq^FF,j,^ma.f"' (4.2) 



where \sq5FF,j,max denotes the solution of eq. (|3.6| ) with the coefficients ( ^.3.1|) -( ^.3.7|) that 
has the maximal magnitude in region Ri. 

As with the cylindrical strips, we can use our new exact solution for the Ly = 5 open 
square strip to study the approach of W to the limit for the infinite 2D square lattice, 
extending ||2^. In Table ^ we list various values of W{sq{Ly x oo, FBCy, BCx),q) (which. 



for this range of q, are independent of BCx), denoted as W{sq{Ly),F,q), together with 
Monte Carlo measurements of W for the full 2D square lattice, W{sq, q) from and 
the g = 3 value W{sq,3) = (4/3)^/2 from §. We also list the ratio Rw{sq{Ly), FBCy,q) 
defined in ( p.l5|) . In ||2^ it was proved that for FBCy the approach of W to the Ly = oo 
limit is monotonic. One sees from Table ^ that for Ly = 5 and moderate values of q, say 
5 or 6, the agreement of W{sq{Ly),q) for the infinite-length, finite-width strips with the 
respective values W{sq, q) for the infinite square lattice is very good, accurate to a few per 
cent, although the approach is somewhat slower for open strips than for cylindrical strips. 
This is understandable since the condition of periodic boundary conditions in the transverse 
direction minimizes finite-size effects in this direction. 
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Table 3: Values of W{sq{Ly x oo, FBCy, BCx,q) (for any BCx), denoted W{sq{Ly), F,q) for short, with 
W{sq,q) = W(sq{po x (x>),q) for 3 < g < 10. For each value of g, the quantities in the upper line are 
identified at the top and the quantities in the lower line are the values of Rw{sq{Ly), FBCy, q). 



Q 


W{sq{l),F,q) 


W{sq{2),F,q) 


W{sq{3),F,q) 


W{sq{4),F,q) 


W{sq{5),F,q) 


W{sq,q) 




3 


2 


1.73205 


1.65846 


1.624945 


1.60597 


1.53960.. 






1.299 


1.125 


1.077 


1.055 


1.013 


1 




4 


3 


2.64575 


2.53800 


2.48590 


2.45517 


2.3370(7) 






1.284 


1.132 


1.086 


1.064 


1.051 


1 




5 


4 


3.60555 


3.48304 


3.42336 


3.38805 


3.2510(10) 






1.230 


1.109 


1.071 


1.053 


1.042 


1 




6 


5 


4.58258 


4.45136 


4.38717 


4.34910 


4.2003(12) 






1.190 


1.091 


1.060 


1.0445 


1.035 


1 




7 


6 


5.56776 


5.43073 


5.36348 


5.32353 


5.1669(15) 






1.161 


1.078 


1.051 


1.038 


1.030 


1 




8 


7 


6.55744 


6.41623 


6.34677 


6.30545 


6.1431(20) 






1.1395 


1.067 


1.0445 


1.033 


1.026 


1 




9 


8 


7.54983 


7.40548 


7.33434 


7.29199 


7.1254(22) 






1.123 


1.060 


1.039 


1.029 


1.023 


1 




10 


9 


8.54400 


8.39720 


8.324745 


8.28157 


8.1122(25) 






1.109 


1.053 


1.035 


1.026 


1.021 


1 
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5 Comparative Discussion on P and B 



In this section we give a general discussion of some properties of (i) the chromatic polynomials 
for cyclic lattice strips with both arbitrarily great length and arbitrarily great width, and 
(ii) the loci B for the infinite-length limit of strips of the square and triangular lattice with 
various boundary conditions. This discussion incorporates the exact solutions given in the 
present work and also in our previous papers. 

1. From our exact solutions for cyclic and Mobius strips of the square and triangular 
lattices, we draw the following inference: for these lattice strips, with arbitrary Ly 
(independent of L^.), the \Gs,j ^q. ( |1.4| ) with the highest-degree c^'^^ namely c^^^-* 
H (see eq. (^)), is 

= (-1)^" (5.1) 

2. A second inference concerns the set of terms Xcsd ^'-'^ cyclic strips of the square 
and triangular lattice with coefficients cg^j = c'-^^^^\ There are Ly of these terms 
^Gs,j il- Let Ag,,l„j = (-l)^''AG„Ly,i for = sq, tri. Then the Xsq,Ly/s and hence 
the \sq,Ly/s with coefficients Cgq^Lyj = c^^y~^^ can be calculated as follows. Denote the 
equation whose solution is Xsq,Ly,j as f{sq, Ly,^). Thus, f{sq, 1, = ^ + (g — 1) and 

f{sq,2,0 = fisq,l,m + q-^) ■ (5-2) 
The Xsg^Ly,/^ for higher values of Ly are then given by 

f{sq,Ly,0 = f{sq,Ly-l,0{^ + q-3)-f{sq,Ly-2,0 for > 3 . (5.3) 

We find that in the chromatic polynomial for the cyclic strip of the square lattice, of 
the Asgj's with coefficient Csqj = c^^i'~^\ (i) one is Xsqj = (— 1)'^''(1 — q); (ii) if Ly is 
even, then another is 3 — g; (iii) if = mod 3, two others are (—1)^^(2 — q) and 
(— 1)'^«(4 — g); (iv) if Lj^ = mod 6, then two others are 3 ± \/3 — q. (This is not an 
exhaustive list of special factors.) 

For cyclic strips of the triangular lattice, denote the equation whose solution is Xtri,Ly,j 
as f{tri, Ly,^). Thus, f(tri, 1, ^) = ,^ + (g — 1) and 

fitrt, 2, = fitrt, 1, e)(e + g - 3) - (e - 1) (5.4) 

The Xtri,LyS^ fo^ higher values of Ly are then given by 

/(trz,L„0=/(tr^,L,-l,0(e + g-3)-e/(tr^,L,-2,0 for > 3 (5.5) 
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The equations defining the Xtri,Ly,i^ involve progressively higher degrees in ^. 



It was shown in [^] that the Ag^j's are the same for the cyclic and Mobius strips 
of a given lattice with width Ly. Therefore, for each width Ly, the A^^j's identified 
above also occur in the respective Mobius strips of the square and triangular lattices, 
although they do not, in general, have the same coefficients cq^j- 

These inferences are important because they show how one can reduce the problem of 
calculating the Ag^j's for larger-width strips graphs of cyclic and Mobius type from 
those for lower widths without recourse to the usual iterative application of the deletion- 
contraction or coloring matrix methods. That is, after having used these latter methods 
to obtain the chromatic polynomials for the first few values of Ly, the rest can be ob- 
tained purely algebraically, without further direct analysis of the graphs involved. Work 
on constructing the recursive formulas for the other Ag^j's is currently in progress. 

3. For the infinite- length limit of a given strip graph Gs, the dominant XgsJ region Ri 
is independent of the longitudinal boundary condition, and its coefficient is c^^^ = 1 
H. In particular, this Ag.j is the same for {FBCy, {T)PBC^) and {FBCy,FBC^) 
boundary conditions. When this XgsJ is the root of an algebraic equation of degree 
higher than linear, then, for the theorem on symmetric functions of roots of algebraic 
equations to apply and guarantee the polynomial nature of P{Gs,q), it is necessary 
and sufficient that all of the other roots of this equation enter with the same coefficient 



[ p9|] . Hence, the analysis given in for cyclic strips of type Gs that determines 
the number of Ag^j's with a specified cqsJ = c^^^ places, for d = 0, an upper bound 
on the number A^a of Ag^j's that occur in the strip of type Gg with {FBGx, FBGy) 
boundary conditions. In particular, for Ly from 1 through 8, this number np{Ly,0) 
takes the values 1,1,2,4,9,21,51,127. For the square lattice, one finds, for Ly from 1 
through the current results presented here for Ly = 5, the values Nx = 1, 1, 2, 3, 7, as 
listed in Table |l]. For the strips of the triangular lattice, this upper bound is realized 
as an equality: for Ly from 1 to 5, the open strips have Nx = 1, 1, 2, 4, 9. The reason 
that the inequality is realized as an equality for the strips of the triangular lattice is 
a consequence of the different behavior of the coefficients of the square and triangular 



lattice strips in the Mobius case . 



4. For all of the strips of the square lattice containing global circuits that we have studied, 
the locus B encloses regions of the q plane including certain intervals on the real axis 
and passes through q = and g = 2 as well as other possible points, depending on 
the family. Note that the presence of global circuits is a sufficient, but not necessary. 
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condition for B to enclose regions, as was shown in ||2^ (see Fig. 4 of that work). Our 
present results for the square lattice are in accord with, and strengthen the evidence 



for, the inference (conjecture) ||3^, Q that 



{q = 0, 2} for sqiLy, FBCy, {T)PBC^) \f Ly > 1 

and sq{Ly, PBCy, {T)PBC^) \/ Ly > 3 . (5.6) 

(For the upper line of this equation, note that the Ly = 1 graphs with {FBCy, TPBCx) 
and [FBCy, PBCx) boundary conditions are identical.) 

5. The crossing of B at the point q = 2 for the (infinite- length limit of) strips with global 
circuits nicely signals the existence of a zero-temperature critical point in the Ising 
antiferromagnet (equivalent to the Ising ferromagnet on bipartite graphs). This has 
been discussed in |33| in the context of exact solutions for finite-temperature Potts 
model partition functions on the Ly = 2 cyclic and Mobius strips of the square lattice. 
In contrast, this connection is not, in general, present for strips with free longitudinal 
boundary conditions since B does not, in general, pass through q = 2. (Of the strips 
with FBCx that we have studied so far, such a crossing at g = 2 was only found for 
the Ly = 3 {FBCy, FBCx) case, as one can see from Table |1|.) Furthermore, for the 
strips without global circuits, there is no indication of any motion of the respective 
loci B toward q = 2 as Ly increases. 

6. Our exact solutions show that in the limit as L^, — oo, the respective loci B for the 
W functions for the infinite strips of the square lattice with the fixed values of Ly 
considered and with (i) periodic or twisted periodic longitudinal boundary conditions 
and (ii) free longitudinal boundary conditions differ; in particular, the loci B for cases 
with (i) pass through q = 2, whereas the loci for cases with (ii) do not. This dependence 
of B on the boundary conditions means that an n — oo limit does not exist in a 
manner independent of these boundary conditions. If one fixes g = 2 at the outset, 
i.e. considers the Ising antiferromagnet on the square-lattice strips (or if one fixes q to 
the trivial value g = 1 at the outset) and then calculates W, there are no pathologies; 
these arise when one considers nonintegral real q in the range < g < 3. This was 
already discussed in the more general context of the full temperature-dependent free 



energy for the Potts antiferromagnet in |33|, together with other pathologies such as a 



negative partition function (lack of Gibbs measure), noted earlier in ||2^, and negative 
specific heat. In general, the the conclusion is that a nonpathological n oo limit 
of the antiferromagnetic Potts model fails to exist at sufficiently low temperature and 
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sufficiently small real nonintegral positive q on strips of the square lattice. Since these 
strips are of fixed width, the oo limit may be considered to be effectively quasi- 

one-dimensional; in contrast, a true two-dimensional thermodynamic limit would be 
Lj. — i> oo, Ly oo, with the ratio Ly/L^ a finite nonzero constant in this limit. 
However, as is clear from the random cluster representation of the Potts model, the 
problem of a negative partition function (lack of Gibbs measure) for sufficiently small 
positive real nonintegral q is present for both infinite-length, finite width strips and 



for the above two- or higher- dimensional infinite volume limit [^, Our exact 

results for infinite-length strips of various widths and our inference above that in the 
Ly ^ oo limit, the loci B and W functions obtained with periodic (or twisted periodic) 
versus free longitudinal boundary conditions would differ is in connected with the 



other pathologies noted above. From the analysis in |^5[, we also conclude that a 
nonpathological ^ oo limit for the antiferromagnetic Potts model fails to exist 
at sufficiently low temperature and sufficiently small positive nonintegral q on strips 
of the triangular lattice and a nonpathological thermodynamic limit fails to exist at 
sufficiently low temperature for nonintegral < g < 4 for the full triangular lattice. 
One could infer a generalization of this for other lattices also: a thermodynamic limit 
would fail to exist for the Potts antiferromagnet at sufficiently low temperature for 
positive nonintegral q in the range from to qc for the given 2D lattice, e.g., gc = 3 
for the kagome lattice. Our exact solutions are consistent with the understanding that 
the point qc for the infinite 2D (or higher-dimensional) lattice is independent of the 
boundary conditions used to define this infinite lattice. 

7. For cyclic strips, we note a correlation between the coefficient cg^j of the respective 
dominant A^^ j's in regions that include intervals of the real axis. Before, it was 
shown that the CGsJ of the dominant XgsJ region Ri including the real intervals 
q > qci{G}) and g < is c^^^ = 1, where the c^'^^ were given in eqs. (|2.8| ), (|2.9|) . 



We observe further that the cg^j that multiplies the dominant XgsJ iii the region 
containing the intervals < g < 2 is c^^\ For the cyclic Ly = 3 and Ly = A strips, 
there is also another region containing an interval 2 < g < g^ on the real axis, where 
g^ ~ 2.34 and 2.49 for Ly = 3,4; in this region, we find that the Cq^j multiplying the 
dominant XgsJ is c^^-*. 

8. Our new results on cylindrical and open strips with Ly = 5 con&m and extend various 



features that had been discussed earlier p3|, |25[: for these values of Ly, B forms 



arcs, and as Ly increases, these arcs elongate and move closer together, with the arc 
endpoints nearest to the origin moving toward this point. This is consistent with the 
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inference that in the Ly oo hmit, the arcs would close to form a closed boundary 
that contained q = and q = qdsq) = 3. One sees this general trend in the Ly = Q 
cylindrical strip (Fig. |]). However, in contrast with the strip graphs containing global 
circuits, for which the loci B contained a region-enclosing boundary passing through 
g = for any Ly, this feature is evidently only approached in the limit as — >^ oo 
for the strips that do not contain global circuits. The earlier calculations of cylindrical 
strips of the triangular lattice showed an example of a strip, namely the Ly = A case, 
where B contains arcs and an self-conjugate oval on the real axis [^, but for the 



cylindrical strips of the square lattice that we have investigated so far, we have not yet 
encountered such an oval. 

For the L^ — > oo limit of all of the strips of the square lattice containing global 
circuits, a qc is defined, and our results for the cyclic and Mobius strips with widths 
from Ly = 1 through Ly = 4 indicate that qc is a non-decreasing function of Ly in 
these cases. The same behavior was found for the strips of the triangular lattice with 
Ly = 2 1^ (and subsequently also Ly = 3,4 ||3^). This motivated the inference 



(conjecture) that qc is a non-decreasing function of Ly for strips of regular lattices with 



(FBCy, {T)PBCx) boundary conditions and our present results strengthen the 



support for this inference. Given that, as — > oo, gc reaches a limit, which is the qc 
for the 2D lattice of the specified type (square, triangular, etc.), this inference leads to 
the following inequahty: 

qc{A, Ly X oo, BCy, {T)PBCx) < q^A) . (5.7) 

Our exact solutions show that this inequality can be saturated. For example, qc = 3 
for the Ly = 3 torus and Klein bottle strip of the square lattice [0, which is equal to 
the qc value for the infinite 2D square lattice [^. In contrast, for (the L^. — ^ oo limit 
of) strips without global circuits, the locus B does not necessarily cross the real axis, 
and hence there is not necessarily any qc defined, as was shown in [Q. Furthermore, 
in these cases, even if B does cross the real axis, so that a qc is defined, the value of qc 
is not a non-decreasing function of Ly. This is shown by our calculations of B for the 
Ly = 4, Ly = 5, and Ly = 6 strips of the triangular lattice with cylindrical boundary 



conditions in [^]; for these we get = 4 for Ly = A but qc = 3.28 for Ly = 5 and 
qc = 3.25 for Ly = 6. Similarly, for the Ly = 5 and Ly = 6 cylindrical strips of the 
square lattice we get qc = 2.69 and q = 2.61, respectively. 

10. A generalized conjecture would be to consider a slab of a ci-dimensional lattice A of 
size Li X L2 X ... X Ld, and let d — 1 of the lengths of this slab go to infinity, holding 
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one length, which can be chosen without loss of generality to be L^, fixed and finite, 
and to define W via ( tL.2|) as 

WiA,LdXoo''-\BCi,...,BCd,q)= lim P(A, LiX...xL,, .., g)^/" 

Lj^oo, j=l,...,d—l 

(5.8) 

For each of these W functions, one would consider the corresponding continuous singu- 
lar locus B and its Qc, for choices of the BCj and where this point exists. We display 
the dependence of qc on these inputs by writing it as qd-^d, Ld x oo'^^^, BCi, BCd)- 
Next, we define a W function for the rf- dimensional lattice as 

Ty(A,,5Ci,...,SC,,g) = lim P(A, Li x ... x L,, 5Ci, 5C,, g)^/" (5.9) 

Lj—KX), j=l,...,d 

and a corresponding singular locus and qc{A.d). As indicated in the notation, one 
expects that this qc would be independent of the BCj, j = 1, d just as is the case for 
the exactly known q^ values for certain 2D lattices. Then we conjecture the inequality 

q^Ad, Ld X oo'^-i, ECi, BCd) < qd^d) (5.10) 

Similarly, a generalization of our inference that qc is a non-decreasing function of 
Ly for the strips with {FBCy, (T)PBCx) would be the conjecture that qc{Ad,Ld x 
fy^d-i^ (T)P_BCi, (T)PBCd-i, FBCd) is a non-decreasing function of Ld- Our exact 
solutions for strips with {PBCy, FBC^) boundary conditions show that if one uses pe- 
riodic rather than free boundary conditions in the direction in which the slab is finite, 
then the resultant qc is not, in general, a non-decreasing function of Ld- 

11. Our exact solutions for the Ly = 4 cyclic and Mobius strips of the square lattice yield 
a singular locus B that has support for Re{q) < 0. In comparison (see Table |l]), this 
was also true for the same type of strip with Ly = 3, while for Ly = 1,2, B only had 
support for Re{q) > 0, and the only point on B with Re{q) = was q = itself. This 
shows that for a given type of strip, increasing Ly can shift the left-most chromatic 
zeros and, in the — > cxd limit, the left-most portion of the locus i3,into the Re{q) < 
half-plane. The same type of behavior was found for the cyclic and Mobius strips of 
the triangular lattice; for Ly = 2 and Ly = 3, B and chromatic zeros had support only 
for Re{q) > 0, while for Ly = 4, this support extended into the Re{q) < region. 
In p3| it was conjectured that global circuits were a necessary condition for lattice 
strips to have chromatic zeros and, in the limit — > oo, a locus B with support 
for Re{q) < 0. However, this conjecture was ruled out by our exact solutions for 



23 



chromatic polynomials, W, and B for homeomorphic expansions]^ of lattice strips with 
{FBCy, FBCx) boundary conditions in |2^, as also by the results for lattice strips 
with {PBCy, FBCx) in [l^. The homeomorphic expansions in have the effect of 



increasing the girth of these strip graphs, and it was found that for a given type of 
open strip graph, increasing the degree of homeomorphic expansion and hence the girth 
shifts the left-most chromatic zeros and, in the limit — oo, the left-most portion of 
B, farther to the left. This is thus a different way of getting chromatic zeros and part 
of B to have support for Re{q) < than in the present case of cyclic strips, where this 
result is obtained as a consequence of increasing the width of the strip while the girth 
remains constant. We remark that for all of these families of graphs, the magnitudes of 
the chromatic zeros and points q on B are bounded. Yet another way to get chromatic 
zeros and B with negative real parts involves families with unbounded chromatic zeros 
and loci B |3^; indeed, in we constructed families where these zeros and loci 



B had arbitrarily large negative Re{q) 

12. There have been a number of theorems proved concerning real chromatic zeros. An 
elementary result is that no chromatic zeros can lie on the negative real axis q < 0, 
since a chromatic polynomial has alternating coefficients. It has also been proved that 
there are no chromatic zeros in the intervals < g < 1, and 1 < g < 32/27 [|^. The 



bound of 32/27 in |^0[ has been shown to be sharp; i.e., for any e > 0, there exists a 



graph with a chromatic zero at q = 32/27 -|-e Based on our studies of strips of the 
square (and triangular) lattices with all of the various boundary conditions considered, 
we make the following observation: for such strips, we have not found any chromatic 
zeros, except for the zero at g = 1, in the interior of the disk |g — 1| = 1. This motivates 
the conjecture that for these strips, there are no chromatic zeros with |g — 1| < 1 except 
for the zero at g = 1. Assuming that this conjecture is valid, the bound would be a 
sharp bound, since the circuit graph with n vertices, C„, has chromatic zeros lying 
precisely on the circle | g — 1 1 = 1 and at g = 1 [0 . 

6 Values of W for Low Integral Values of q 



In previous works |26|, ^ and sections of the present paper we have discussed values of 
W for various infinite-length, finite-width lattice strips. For infinite-length limits of strips 
with global circuits, where the region(s) of the positive real axis in the interval < g < gc 

^We recall two definitions from graph theory: (i) a homeomorphic expansion of a graph is obtained by 
inserting one or more degree-2 vertices on edge(s) of the graph; (ii) the girth of a graph is the number of 
edges or vertices in a minimum- distance circuit. 



24 



are not analytically connected with the region Ri including q > qc (and g < 0), the ground 
state degeneracy per site, W, has a qualitatively different behavior than for integer or real 
q > qc- A comparative discussion of this was given in with the results available at that 
time, and it is worthwhile to use our new exact solutions to study this behavior further. In 
particular, it is of interest to inquire what the values of the W functions are for the infinite- 
length limits of various strips of the square lattice at the points g = 0, 1, and 2. Our exact 
analytic expressions yield the numerical values listed in Table ^. The notation follows that 
in Table |I[ As was noted in 0, in general, for regions other than Ri, it is only possible 
to determine \W\ unambiguously. Hence, for uniformity, we list \W\ for all of the strips, 
including those with only a region Ri. For comparison, we also include values of \W\ at 
g = 0, 1, 2 and 3 for infinite- length strips of the triangular lattice in Table |^. In addition, for 
families where, in the Lr,. —>■ oo limit, there exists a qc, we include the respective values of 
W at qc- For the smallest widths, the \W\ values are relatively simple analytic expressions, 
e.g., for the square strips, \W\ = S^/^ fo^ {FBCy,BC^), Ly = 2, q = 0; \W\ = 13^/^ for 
{PBCy, BCx), Ly = 3, q = 0, and so forth. In the case of the triangular lattice, Ly = oo, 
[PBCy, FBCx), the values of \W\ for g = and g = 4 are from [0; the exact value for 
g = 3 is our analytic evaluation, and the numerical values for g = 1, 2 are our numerical 
evaluations, of an integral representation in Although we list the values in the tables 
only to three significant figures, we note that the g = 1 value, |W^(trz)| ^ 3.1716, is different 
from the g = value \W{tri, 3 x oo, PBCy, BC^)] - 3.1748. 



For these values of g, the noncommutativity of eq. (|1.3| ) occurs [Q. Thus, for any 
connected graph G, and in particular, the lattice strips considered here, the chromatic poly- 
nomial P{G, q) vanishes at g = and g = 1 and hence also the function Wnq defined via 
the order of limits on the right-hand side of eq. (|1.3| ) vanishes. In contrast, in general, Wqn 
defined by the limits on the left-hand side of ( |1.3|) is nonzero. For cyclic strips of the square 
lattice of length L^, at g = 2, the chromatic polynomial P is equal to 2 if L^ is even but 
if Lx is odd, so that at g = 2, no Wnq is defined, since the limit on the right-hand side 
of ( |1.3|) does not exist; however, Wqn is well-defined and, in general, nonzero. Analogous 
comments apply for strips of the triangular lattice: at g = 2, the chromatic polynomial P 
vanishes identically, so Wnq = 0, but Wqn is, in general, nonzero. For cyclic strips of the 
triangular lattice, at g = 3, then P = 3! if = mod 3, and P = if L^. = 1 or 2 mod 3; 
hence, no Wnq is defined, since the limit on the right-hand side of ( |1.3D does not exist, but 
Wqn is well-defined and, in general, nonzero. As with the other results given in this paper, 
the values of W given in Tables H and || follow the definition W = Wgn- 

Some general comments follow: 
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Table 4: Values of W{sq,Ly x oo, BCy, BCx,q) for fow integral q and for respective qc, if such a point 
exists, where BCy and BC^ denote the transverse and longitudinal boundary conditions. 



Ly 


BCy 


BC, 


\Wg=0\ 


\W,=i\ 




Qc 




1 


F 


F 


1 





1 


n 




2 


F 


F 


1.73 


1 


1 


n 





3 


F 


F 


2.06 


1.44 


1 


2 


1 


4 


F 


F 


2.24 


1.68 


1.25 


2.30 


1.14 


5 


F 


F 


2.36 


1.82 


1.39 


2.43 


1.22 


1 


F 


P 


1 


1 


1 


2 


1 


2 


F 


(T)P 


1.73 


1.41 


1 


2 


1 


3 


F 


(T)P 


2.06 


1.66 


1.26 


2.34 


1.18 


4 


F 


(T)P 


2.24 


1.81 


1.41 


2.49 


1.27 


3 


P 


F 


2.35 


1.59 


1 


n 




4 


P 


F 


2.58 


1.89 


1.32 


2.35 


1.16 


5 


P 


F 


2.68 


2.05 


1.51 


2.69 


1.15 


6 


P 


F 


2.73 


2.14 


1.62 


2.61 


1.39 


3 


P 


(T)P 


2.35 


1.91 


1.44 


3 


1.26 



Table 5: Values of W{tri,Ly x oo, BCy, BCx,q) for low integral q and for respective qc, if such a point 
exists, where BCy and BCx denote the transverse and longitudinal boundary conditions. 



Ly 


BCy 


BCx 






\W,=2\ 




Qc 




2 


F 


F 


2 


1 





1 


n 




3 


F 


F 


2.49 


1.66 


1 


1 


2.57 


0.656 


4 


F 


F 


2.77 


2.02 


1.41 


1 


n 




5 


F 


F 


2.95 


2.25 


1.66 


1 


3 


1 


2 


F 


(T)P 


2 


1.62 


1 


1 


3 


1 


3 


F 


(T)P 


2.49 


1.99 


1.44 


1 


3 


1 


4 


F 


(T)P 


2.77 


2.22 


1.68 


1.21 


3.23 


1.13 


3 


P 


F 


3.17 


2.22 


1.26 


1 


n 




4 


P 


F 


3.44 


2.60 


1.78 





4 


1.19 


5 


P 


F 


3.56 


2.79 


2.05 


1.15 


3.28 


0.772 


6 


P 


F 


3.63 


2.90 


2.21 


1.41 


3.25 


1.11 


oo 


P 


F 


3.77 


3.17 


2.60 


2 


4 


1.46 


3 


P 


(T)P 


3.17 


2.62 


2 


1.71 


3.72 


1.41 
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1. As is evident in Tables ^ and ^ for values of q that are positive but sufficiently small, for 
a given lattice, boundary conditions, and value of Ly studied, \W\ is a non-increasing 
function of q. In contrast, for sufficiently large g, \W\ increases with q. For families 
of graphs that involve global circuits, these two different types of behavior occur, 
respectively, for < g < gc and q > qc- The latter behavior is the one expected for the 
g-state Potts antiferromagnet, since increasing g increases the physical ground state 
entropy. As examples, for the {L^ —>■ oo limit of) circuit graph, W is constant for 
< g < 2, while for the cyclic or Mobius strip of the square lattice Lj^ = 2, it decreases 
as \W\ = |3 — g| in this interval; in both of these cases, gc = 2 and W is real and 
increasing for q> q^.. For the cyclic or Mobius Ly = 3 and Ly = 4 strips of the square 
lattice, \W\ has a different analytic form in the interval < g < 2 and 2 < q < q^ but 
is everywhere decreasing for < g < gc, for the respective values of qc- As an example 
of a strip with no gc, for the open line, Ly = 1, \W\ decreases from 1 to as g increases 
from to 1 and increases for larger g. As another example of a strip with no gc, for 
the Ly = 3 strip of the square lattice with {PBCy, FBCx) boundary conditions, \W\ 
decreases monotonically as g increases from and vanishes at g ~ 2.453; for larger 
values of g, W is real and positive and increases with g. 

2. For the L^ ^ oo limit of strips with free transverse boundary conditions, FBCy and 
any longitudinal boundary conditions BCx-, it was proved that for a fixed physical 
g > gc, is a monotonically decreasing function of Ly |2^. However, as is evident 
from Tables ^ and |^ for sufficiently small positive values of g (smaller than gc for strips 
with a gc), \W\ is a non- decreasing function of Ly. 

3. For the strips that we have studied whose ^ oo limit yields a locus B with a gc, 
|W^(g)| for fixed g G [0, gc] is a non- decreasing function of Ly. 

4. It has been shown that for physical values of g in the g-state Potts antiferromagnet, in 
the La; ^ oo limit of a strip of a given type of lattice A, W^(A, Ly x oo, BCy., BCx-, q) 
is independent of the longitudinal boundary condition BCx However, for small 
positive values of g, \W\ does depend on both BCy and BCx, as is evident from 
Tables ^ and ^. One observes that for the small integral values of g shown in these 
tables, \W{A,Ly x oo, FBCy, F BCx, q)\ < \W{A,Ly x oo, FBCy, {T)PBCx,q)\ and 
\W{K,LyX oo,PBCy,FBCx,q)\ < \W{A, Ly x oo, PBCy, iT)PBCx, q)\. 

5. It was observed in and proved in (section 7 of) p2| that for integer, and, by analytic 
continuation, real, values of g > max{qc) for the square and triangular lattice strips, 
i.e., g > 4, W(tri,q) < W{sq,q). Most of the values of \W\ shown in Tables § and || 
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show the opposite inequahty. Together with various other properties noted above, this 
shows that \W\ behaves quahtatively differently for sufficiently small positive values of 
q than for larger values. 

7 Conclusions 

In conclusion, we have presented exact solutions of the zcro-tcmpcraturc partition function 
(chromatic polynomial P) and the ground state degeneracy per site W (= exponent of the 
ground-state entropy) for the g-state Potts antiferromagnet on strips of the square lattice 
of width Ly vertices and arbitrarily great length vertices. The specific solutions were for 
(a) Ly = 4, {FBCy,PBC^) (cyclic); (b) Ly = 4, {FBCy,TPBC^) (M5bius); (c) Ly = 5,6, 
{PBCy,FBC^) (cylindrical); and (d) Ly = 5, {FBCy,FBC^) (open), where FBC, PBC, 
and TPBC denote free, periodic, and twisted periodic boundary conditions, respectively. 
Some inferences were given for certain terms Ag^j for cyclic and Mobius strip graphs of the 
square and triangular lattice that allow one to calculate them for arbitrarily wide strips (of 
any length). These are important because they show how one can reduce the problem of 
calculating the Ag^j's for these strips of arbitrarily large width from those for lower widths 
without recourse to the usual iterative application of the deletion-contraction or coloring 
matrix methods. A comparative discussion was given of the continuous nonanalytic locus B 
for these strips and numerical results of W were given for a range of values of q. In general, 
our exact solutions give further insight into the properties of the Potts antiferromagnet in 
the setting of infinite-length, finite width systems. 

Acknowledgment: The research of R. S. was supported in part by the NSF grant PHY- 
9722101 and at Brookhaven by the DOE contract DE-AC02-98CH10886.6 

8 Appendix 

8.1 Terms Aj for the Cyclic 4 x m Strip of the Square Lattice 

In this appendix we give the equations for the terms \j for 7 < j < 26. The Xsq4j, j = 7, 8, 9, 
are roots of the cubic equation 

+ + V - 23?" + Alq - 33)f 

^Accordingly, the U.S. government retains a non-exclusive royalty-free license to publish or reproduce the 
published form of this contribution or to allow others to do so for U.S. government purposes. 
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+(2g^ - 23g^ + 116g^ - 329g^ + 553g^ - 517? + 207)^ 

+ {-q^ + 16g^ - 112g^ + 449g^ - llSOg^ + 1829g3 - 1858?^ + 1084g - 279) - . 

(8.1.1) 

The Xj for j — 10, 11, 12 are the roots of another cubic equation 

e - 2{q -2){q- 3)^' + {q" - llq' + 42^^ - 68? + 38)^ 

-{q - l)(g - 3)(-g2 + 7g2 - 15q + 11) = . (8.1.2) 
The Xj for 13 < j < 16 are the roots of the quartic equation 
+ (2q^ - 12q^ + 28q - 23)^^ 

+ {q^ - 13q^ + 73g^ - 224g3 + 396?^ - 381? + 152)^^ 

-{q-l){q- 3){q^ - I2q^ + 62g^ - 179g^ + 304^^ - 288g + 119)^ 

-q^ + 17g^ - 127g^ + 549g^ - 1518g^ + 2790g^ - 3411g3 + 2673^^ - 1215g + 243 = . 

(8.1.3) 

Finally, there are two sets of roots of two degree-5 equations. The set Xj for 17 < j < 21 are 
the roots of the equation 

^' + (q-3){2q'-9q + U)e 

+ {q^ - 17q^ + 119g'^ - 446g^ + 947^^ - 1080g + 515)^^ 

-{q - 3)(2g^ - 32g^ + 224g5 - 883g^ + 2106g^ - 3028^^ + 2428g - 835)^^ 

+ {q - 2){q^ - 22q^ + 213/ - 1193^ + 4267g^ - 10120g^ + 15922g^ - 16013?^ + 9329g - 

+ {q - l){q^ - 17/ + 125/ - 520g^ + 1342/ - 2206g^ + 2261/ - 1325g + 341) (g - 3)^ = 
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(8.1.4) 

The set Xj for 22 < j < 26 are the roots of the equation 

^5 + (_4g2 ^ _ 2g^^4 ^ (^g^4 _ 5g^3 ^ 214g2 - 354g + 219)^^ 

+ (-4g^ + 60g^ - 370g^ + 1198g^ - 2144^^ + 2013g - 773)^^ 

+ {q - 2)(g^ - 20q^ + 162g^ - 693g^ + 1697g^ - 2391g^ + 1805g - 565)^ 

+ (g - l)(g - 3)(g^ - 16g^ + 106g^ - 378g^ + 788g^ - 967^^ + 653g - 189) = . 

(8.1.5) 

8.2 Generating Functions for the Ly = 5,6 Strips of the Square 
Lattice with {PBCy, FBC^) 

For the Ly = 5 strip we calculate a generating function of the form ( p.2| ) with dj) = 2, d_\f = 1 
and, in the notation of eqs. (|3.4|) and ( |3.3D , we find 

bsg5PF,i = -q^ + lOg^ - 46g3 + 124^^ - 198g + 148 (8.2.1) 
bsq5PF,2 = q^- 19g^ + 159g'' - 767g^ + 2339g^ - 4627g^ + 5800^^ - 4212g + 1362 (8.2.2) 
Asq5PF,o = q{q - l)(g - 2)(g^ - 12g^ + 67g^ - 225g^ + 494g^ - 719g^ + 650g - 282) (8.2.3) 

Asqt>PF,i = -q{q - l)(g - 2)(g' - 2g + 2)(g« - 19g^ + 159g'' - 767g^ 

+2339g^ - 4627g^ + 5800g2 - 4212g + 1362) (8.2.4) 

For the Ly = Q strip we calculate a generating function of the form ( |3.2|) with = 5, 
dj\f = 4, with 

bsg6PF,i = -q^ + I2q^ - 68g^ + 234g3 - 524^^ + 727g - 483 (8.2.5) 
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bsq(iPF,2 = 2g^° - Uq^ + 456g^ - 2917g^ + 12710g^ - 39322^5 + 87323g^ - 137193g3 
+14562V - 94100? + 28114 (8.2.6) 

bsq6PF,3 = -q^^ + 33g^^ - 509?^^ + 4872g^^ - 32374?^° + 158152^^^ - 586234?^ 

+1676100g^ - 3715937g^ + 6358772g^ - 8268225g^ + 7921161?^ 

-5284418g^ + 2197026g - 429510 (8.2.7) 

bsq(iPF,4 = + 38g^^ - 681?^^ + 7649g^^ - 60357?^^ + 355400?^^ - 1618550g^^ 
+5828269?^° - 16812727?^ + 39098146?^ - 73327191g^ + 110295876?^ - 131415610?^ 
+121386275g^ - 83893487g^ + 40850378?^ - 12502528g + 1809361 (8.2.8) 

bsqePF,5 = q^^ - 41g^^ + 794g^^ - 96585^^^ + 82760g^^ - 531052g^^ 

+2647330g^3 - 10495556g^2 + 33592560^^^ - 87588439g^° + 186851845g^ - 326185418g^ 
+464098186g^ - 533530852g^ + 488389118g^ - 347889815g^ + 185960167g^ - 70211630g^ 
+16703951g - 1884267 (8.2.9) 

Asq6PF,o = qiq - l)(g^° - iV + 136g^ - 674g^ + 2296g^ - 5640g^ + 10183g^ - 13457g^ 
+12563?^ - 7517g + 2183) (8.2.10) 

Asq6PF,i = -q{q - l)(2g^^ - 5Aq^^ + 695g^^ - 5631g^^ + 31999g^° - 134668?'' 
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+432404?^ - 1075802g^ + 2085064?^ - 3137110?^ + 3615627?^ - 3106751?^ + 1890461?^ 
-733250g + 137516) (8.2.11) 



Asq&PF,i = q{q - l)(g^^ - 38?^^ + 684?^^ - 7756?^^ + 62128?^^ - 373554?^^ 
+1748131?^^ - 6513823g^^ + 19602672g^° - 48032023?^ + 96128905?^ - 156920332g^ 
+207640116g^ - 220043849g^ + 183010634g^ - 115543495?^ + 52290297?^ 
-15182726g + 2135038) (8.2.12) 

Asq%PF,z = q{q - l)(g^^ - 43g^° + 881?^^^ - 11444?^^ + 105796?^^ - 740641g^^ 
+4078480g^^ - 18111664g^^ + 65961019g^^ - 199240735?^^ ^ 502713558g^^ 

-1063474616g^° + 1887470282^^ - 2803761470g^ + 3465937164g^ - 3530769703g^ + 2919336052g^ 
-1914246633g^ + 960052617g^ - 346744827g^ + 80479446g - 9033772) (8.2.13) 

Asq&PF,A = -q{q - l)(g^ - + lOg^ - lOg + 5)(g^^ - Alq^^ + 794g^^ - 9658g^^ 
+82760g^^ - 531052g^^ + 2647330g^^ - 10495556^^^ + 33592560g^^ - 87588439g^° 
+ 186851845g^ - 326185418g^ + 464098186g^ - 533530852^^ + 488389118g^ 
-347889815^^ + 185960167g^ - 70211630^2 + 16703951g - 1884267) (8.2.14) 



8.3 Generating Function for the Ly = 5 Open Strip of the Square 
Lattice 

For this strip we calculate a generating function of the form (|3.2| ) with dv = 7 and dj^ = 6. 



In the notation of eq. (|3.4| ) we find 

bsq{5),i = -q^ + - 40g^ + lOTg^ - 167g + 118 (8.3.1) 

bsg5FF,2 = ^q^ " 63g^ + 458g^ - 2011g^ + 5840g^ - 11477g^ + 14844g2 - 11466g + 4003 

(8.3.2) 

fe.q5FF,3 = -6g" + 136g^° - 1432g^ + 9250g^ - 40749g^ + 128594g^ 
-296624g^ + 499762g^ - 601803g^ + 492117g2 - 245164g + 56113 (8.3.3) 

hsq^FFA = 4:q^^ - 120g^^ + 1685?^^ - 14681g" + 88695g^° - 393187g^ + 1319323g*^ 

-3404712g^ + 6790667g^ - 10414582g^ + 12084263g^ - 10278730g^ 

+6051725g2 - 2204111g + 373840 (8.3.4) 

bsg5FF,5 = "(? - l)iq^^ - 38g^^ + 674g^^ - 7419g^^ + 56807?^^ - 321258g^^ 
+1389731g^° - 4696189g^ + 12540817g* - 26576855g^ + 44582788g^ - 58613690g^ 
+59234653g^ - 44497390g^ + 23436736?^ - 7733009g + 1204091) (8.3.5) 

bsg5FF,6 = "(? - l)^(g^^ " 38g^^ + 682g^^ - 7680g^^ + 60795g^^ - 359135?^^ + 1639962^^^ 
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-5915021g^° + 17065698?'^ - 39623309?^ + 74056302?^ - 110813572?^ 
+131155616g^ - 120231650g^ + 82455281?^ - 39872376?^ 
+ 12141916g- 1753922) 



(8.3.6) 



bsq5FF,7 = (g - l)^(g - 2)^(g^^ - 34g^^ + 538g^^ - 5259g^^ + 35541g^^ - 176036g^° 
+660682g^ - 1914798?^ + 4324155g^ - 7615130g^ + 10381339g^ - 10768339g^ 

+8235159g3 - 4388527?^ + 1459163g - 228580) . (8.3.7) 

Since the Agg^ppj are rather lengthy, they are given in the copy of this paper in the cond-mat 
archive. With the definition Asq^ffj = q{q — l)^sg5FF,j, we have 

As,5FF,o = (D^)^ (8.3.8) 

where D4 = — 3q + 3; 



Asq5FF,i = -4g^^ + 75q^^ - 653q^ + 3478?^ - 12572?^ + 32346?^ 

-60381g^ + 81687?^ - 78370g^ + 50664?^ - 19788g + 3517 (8.3.9) 

Asq5FF,2 = 6g^^ - 154g^^ + 1854?^^ - 13864?^^ + 71883g^° 

-273164g'' + 784036?^ - 1725384?^ + 2923023?^ - 3789945g^ + 3697547g^ 

-2627998g=^ + 1283656?^ - 384667? + 53170 (8.3.10) 

Aq5FF,3 = -4g^^ + 132^^6 - 2056g^^ + 20067g^^ - 137414?^^ 
+700413?^^ - 2750993?^^ + 8502143?^° - 20926266?'' + 41238149?^ 
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-65036748?^ + 8157462^ - 80321984?^ + 60731068?^ - 34014621g^ 
+13280602g2 - 3222200g + 365089 (8.3.11) 



Asg5FF,4: = (g - if{q^^ - 40?^^ + 751?^^ - 8804?^^ + 72287?^^ - 441816?^=^ 
+2084720^^2 - 7769759g^^ + 23199424?^° - 55934061?'' + 109187879?^ 
-172199452g^ + 217795442g^ - 217905554?^ + 168626444g^ - 97343549?^ 
+39444397g2 - 10000178g + 1191823) (8.3.12) 

Asq5FF,5 = (? " lf{q^^ " 40g^^ + 759g^^ - 9082?^^ + 76836?^^ - 488373g^^ 
+2418556g^^ - 9549855?^^ ^ 30509903g^^ - 79556288g^° + 169998158g^ 
-297636195g^ + 425129474g^ - 490934225^^ + 451509111g^ - 323033449g^ 
+ 173271795g^ - 65533774^2 ^ 15574281g - 1747588) (8.3.13) 

Asq5FF,6 = -{q - if{q - 2)2(g^^ - 34g^^ + 538q^^ - 5259q^^ + 35541g^^ 
-176036g^° + 660682g^ - 1914798g^ + 4324155g^ - 7615130^*^ + 10381339g^ 
-10768339g^ + 8235159g2 - 4388527g2 + 1459163g - 228580) (8.3.14) 
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